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Abstract 

The particle emission in relativistic hydro-dynamic model is formulated assuming a 
sharp 3-dimensional space-time freeze-out hypersurface. The boundary conditions cor- 
respond to the energy-momentum and charge conservation between fluid and the gas 
of free particles . The emission and possible feedback of particles are included into the 
self-consistent hydrodynamic scheme. For the time-like parts of the freeze-out hyper- 
surface the obtained results are different from the well known Cooper-Frye formula and 
lead to the shock-like freeze-out. A simple-wave hydrodynamic solution is considered 
in detail to illustrate some aspects of new freeze-out procedure. 
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1. Introduction 



During last three years an essential progress has been achieved in the understanding 
of the freeze-out problem in relativistic hydrodynamics. A generalization of the well known 
Cooper-Frye (CF ) formula |3J] to the case of the time-like freeze-out hypersurfaces (hereafter 
FHS) was suggested in Ref. 0. This formula does not contain negative particle number 
contributions on time-like FHS appeared in the CF procedure from those particles which 
cannot leave the hydrodynamical system during its expansion. The new distribution function 
has the cut-off factor, and we'll call it the cut-off distribution function (CO). Once the 
CO distribution function is defined, the energy-momentum tensor of the gas of free particles is 
known. In contrast to the CF formula the expressions derived in Ref. |2| show the explicit 
dependence on the parameters of the time-like FHS . 

The freeze-out procedure of Ref. [@ has been further developed in a series of publications 
P, ||, ||, H [7]. However, many important questions have not yet been studied. The aim of 
the present paper is to clarify the connection between new freeze-out procedure of Ref. [[| 
and the fluid hydrodynamics. We show that the system evolution is consistently split into 
the differential equations for the fluid motion and particle emission from FHS : smooth 
CF freeze-out on its space-like parts and the shock-like freeze-out with CO formula on time- 
like ones. The conservation laws and consistency of the full system of equations are shown. 

High energy nucleus-nucleus collisions presumably give us the possibility to study strongly 
interacting matter in the framework of relativistic hydrodynamics (see |8|, |9], [T(| and refer- 
ences therein). It is evident that the hydrodynamic description of the fluid created in A+A 
collisions could be meaningful in the finite space-time region only. The relativistic hydrody- 
namics has to be combined self-consistently with the free stream of final particles towards 
the detectors at the latest stage of A+A reaction. 

In Section 2 several transparent examples are presented to show the difference between 
the results obtained with the traditional CF formula and those with CO distribution. The 
form of the distribution function of the gas of free particles is discussed accounting for both 
the emission and (possible) feedback of the particles from the various time-like FHS . 

In Section 3 we study all possible boundary conditions between the perfect fluid and the 
gas of free particles . The fluid equations of motion are derived there. The global energy- 
momentum and charge conservation holds for the whole system consisting of the fluid and 
the gas of free particles . 

The equations for the time-like parts of the FHS are studied in Section 4. The general 
scheme of their solution together with the hydrodynamical equations for the fluid is dis- 
cussed and applied to hydrodynamical motions in 1+1 dimensions. Thermodynamical and 
mechanical stability of the freeze-out shock is discussed there. The freeze-out problem of the 
simple wave for the massless gas of free particles is solved analytically. Momentum spectra 
for the freeze-out of the simple wave are found and compared with those obtained from the 
Cooper-Frye formula. 

In Appendix A we present the systematic analysis of the angular and momentum integra- 
tions for the gas of free particles . Appendix B gives the complete comparison of the Eckart 
and the Landau-Lifshitz definitions of the hydrodynamic velocities for the massless gas of 
free particles. The spatial anisotropy of the energy-momentum tensor in the Landau-Lifshitz 
frame is found. 
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A short resume of the above results and conclusions are given in Section 5. 
2. Particle Spectra on the Time-like Hypersurfaces 



Originally the idea of freeze-out was postulated by Landau [TT|. It states, that if the 
mean free path of the particles of the rarefied gas becomes comparable with the size of the 
whole system, one can assume that the collisions between particles cease to exist on some 
hypersurface and the particles freely move to the detector. Here we shall not consider the 
justification of this assumption, but we shall show that it leads to rather strong restrictions 
on the possible freeze-out procedure. The first, crucial step for showing this is to obtain the 
correct formula for the particle spectra emitted from both space-like and time-like freeze-out 
hypersurfaces. 



2.1. Short Derivation of the Cut-off Formula 



The correct formula for the particle spectra emitted from an arbitrary hypersurface was 
obtained in [g] by the method developed earlier fl2[]. However, we would like to rederive it 
here in a less rigorous, but more transparent way. In what follows we shall use the Landau- 
Lifshitz convention for the metric tensor and 4-intervals. Thus ds 2 > (ds 2 < 0) is a 
time-like (space-like) interval. In contrast to Ref. [|J we shall specify the hypersurfaces by 
their intervals. 

First, let us consider a small element dS± of the boundary between gas and vacuum 
and let us count the number of particles dN that will cross this element during the time dt, 
assuming that the mentioned boundary is at rest (this assumption is very essential because it 
is valid only for the time-like hypersurfaces ). Let 4> ( x > t, p) be a distribution function which 
describes the probability to find a particle with 3-momentum p and energy po = y/m 2 + p 2 
at the space-time coordinate (x, t). Let n denote a unit 3- vector of the external normal to 
the surface element dS±. Then the desired number of particles is given by the expression 



dN = 0(x,t,p)d 3 pv ■ndtdS±e(v ■ n) . (2.1) 

The presence of the B-function in the above formula is evident: it ensures that only particles 
with positive projection of the velocity on the external normal can cross the surface (cosine of 
the angle between the vectors v and n must be positive). Particles with negative projection 
of the velocity on the normal never cross the surface and should not be taken into account. 
Next we take into account the shift of the surface <ix with the time as follows: 



dN = (f) (x, t, p) d 3 p [vdt - dx] ■ n dS±Q Qydt - dx] ■ n) . ( 2.2) 

The second term in the last equation describes the change of the flux through the element 
dS± due to the displacement <ix of the surface element in time dt (it is still not the most 
general choice of the FHS ). Then it can be written in the fully covariant form, 

Po 33 = 4> ( x > t, P) p v da u 9 (p p da p ) , (2.3) 
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where da^ is the external normal 4- vector to the freeze-out hypersurface. This relation differs 
from the Cooper-Frye one by the presence of the G-function which plays an important role 
for the time-like FHS , as we saw it during the derivation. On the other hand it completely 
agrees with the CF formula M for the space-like FHS . In the latter case the 0-function is 
equal to unity because of the inequality p u da u > 0. Thus, Eq. (2.3) is the CF formula [j], 
but without negative particle numbers that appear for the time-like FHS ! 



2.2. Consideration of the Sun-like Object 

As a simple application of the derived formula let us consider the Sun-like object^ (here- 
after "the Sun," for simplicity), i.e., spherical object with the constant radius R Q , ^jf- = 0, 
with the phase space distribution function of photons / 7 (fr), which is isotropic in the mo- 
mentum space and is spherically symmetric in the coordinate one. Here po is the photon 
energy and T = T(|R|) is the Sun temperature on its radius |R|. 

Then one can easily find the total number of the emitted photons by the CF formula. 
With the help of the following expression for the normal 4-vector da^ 

dcTf, = 

one can trivially perform the integration over momenta in spherical coordinates. Indeed, the 
number of emitted photons per unit time and per unit area is given by the CF formula 

£ = f*» f - (t) m = 2 */*° P « ^ (f ) l'^Mco S (e p) , ( 2.4, 

where we decouple the angular and energy integrations because of the spherical symmetry, 
choose the angle 6 P between Z-axis and the vector of 3-momentum as one of the spherical 
angles, and perform a trivial integration over the other spherical angle. 

Then one sees that the integration over all possible values of the angle 8 P gives zero: 




dN° F r+i , . 

dx x = . ( 2.5) 

l 
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dtdS_ 

The calculation of the emitted energy can be done similarly, and again one gets zero for the 
energy emitted by the Sun. In other words, according to the Cooper-Frye formula all stars 
are invisible! 

In contrast to this, our formula (or Eq. (2.3) ) gives the correct answer because in this 
case the angular integral contains the O-function which ensures that only outgoing photons 
are taken into account: 

dN^ ut ~°ff r % r +1 



dtdS i 



£ dOpSm^p) cos(# p )0 (cos(0 p )) ~ J dx xG(x) ^ . ( 2.6) 



1 This nice example was first suggested to one of us, K.A.B., by D.H. Rischke 
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This is one of the best examples showing an inapplicability of the CF -formula for time-like 
hypersurfaces. 



2.3. Other Formulae for the Invariant Spectrum 



In this subsection we would like to comment our generalization of the CF result for the 
case of time-like hypersurfaces. An attempt to modify the CF formula for time-like FHS was 
made by Sinyukov in Ref. |T3|] . This formulation was recently used in Ref. 



In the derivation presented in Ref. [|T^], the decay of a small element of the gas of free 
particles was considered inside the box with the real wall. This is a suitable treatment for 
the space-like FHS [fL2"|. It appears to be, however, not applicable for the time-like one. 



The treatment of the box wall in Ref. |L2] was just an auxiliary trick in order to simplify 
the derivation. It means an introduction of the identical element of the gas behind the 
wall (see Figs. 1-3). Reflection of the particle trajectories from the wall means nothing 
but only an extension of their trajectories behind the wall and emission of those particles 
with corresponding momenta from the symmetric element behind the wall. Counting all 
particles which cross the time-like FHS as it was done in work [[|, one can easily see where 
the reflection from the wall gives a wrong contribution. 

Let us briefly show this. Suppose that in the rest frame of the small element Ax the 
phase space distribution function (x, t, p) is given. Following Ref. [0, we call this frame the 
rest frame of the gas before the decay, or, in short, the reference frame of the gas (hereafter 
RFG). It is convenient to perform the whole consideration in a two-dimensional case for 
the sake of simplicity. 

The first contribution corresponds to those particles with negative momenta that leave 
the element Ax without a reflection from the wall BD in the meaning of the paper (see 
Fig.l) 

' /Al <f)(x,t,p)Axdp x e(-p x ) . (2.7) 



d 2 p±dS± 

where p± is the transverse momentum of the particle. 

Then, let us count the number of particles with negative momenta that are reflected 
from the wall BD. This is equivalent to the consideration of the particles from the element 
BC = Ax (see Fig. 2.). For example, the real trajectory of the particle emitted in the point 
A is AG, but its reflection from the wall BD is EF. This corresponds to the particle emitted 
in the point C with the passage CF. One can easily find out that this procedure gives the 
correct result for the low values of momenta, namely — ^§ = — < ^o, otherwise it is 
wrong. 

The correct contribution is given by the particles with negative momenta from the element 
BH = -fAt (Fig. 3.) 



dN n x 

-<f) (x, t, p) P — At dp x G(-p x ) . ( 2.8) 



d 2 p±dS± ' ' po 

Since the size of the element BH = —'^At depends upon the momentum p x , then for 
the high negative momenta — 1 < ^ < — ^§ = — the contribution from the element 
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Ref. 



— poAt — Ax (Fig. 3.) cannot be obtained by the reflection procedure used in 
13). In fact, there is no 'wall' in the considered physical problem. 



We also have to mention that the energy-momentum tensor derived with the distribution 



function of Ref. |n| is not symmetric, and hence there are difficulties with the conservation of 
the angular momentum. This fact makes it impossible to diagonalize the energy-momentum 
tensor and, therefore, creates a problem in defining the hydrodynamical velocity according 
to the Landau-Lifshitz prescription. 

It is also evident that for the space-like FHS such a reflection procedure gives the correct 



PL <- AB _ Ax 
p° ^ BD At 

The last contribution comes from the particles with positive momenta p x 



result because the inequality — 1 < 



is not fulfilled anymore. 



from the ele- 
ment Ax — -^-At. We should stress, however, that those particles will cross the freeze-out 
hypersurface BD if and only if their velocity ^ is smaller than the derivative of the 
FHS in the t — x plane. Thus, the third term reads as follows (see Fig. 4): 



dN* 



(x, t, p) 



Ax 



At 



dp x Q{p x )Q 



P 



2.9) 



d 2 p±dS± 

In a previous paper || it was shown that Eqs. (2.7) - (2.9) can be reduced after some 
algebra to Eq. (2.3) which we derived at the beginning of this section. 



Another example for the momentum distribution function can be found in Ref. |L5[ and 
it reads as 



dN ± . 
d 6 p 



p^u^ u v da v 



2.10) 



where u v is the 4-vector of the hydrodynamical velocity. It is clear that this formula gives a 
correct result in the absence of the hydrodynamical motion only. 

Thus, we completed the correct generalization of the CF -formula for the case of arbitrary 
FHS . 



2.4. Interpretation of the Cut-off formula 

Here we would like to give an additional interpretation of the CO -formula for the invari- 
ant spectrum. The presence of the G-function in the momentum distribution, which cuts 
off the unphysical contributions of negative particle number, is also important in the coordi- 
nate space. This is schematically shown for the 2-dimensional case (when the hypersurfaces 
become just curves) in Figs. 5 and 6 for convex and concave hypersurfaces, respectively. 

Situation with the convex freeze-out hypersurface (see Fig. 5) is evident since no particles 
can reenter the fluid. Then, for the fixed velocity v x = ^ of the particles measured by a 
detector, the cut-off function gives the integration limits in coordinate space which are 
obtained from the condition p^da^ = 0, or 

_ dx*(t) p x 
v<?= , = — , 2.11 
dt p [) 
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where the freeze-out hypersurface is given by the equation x = x*(t). The above equation 
gives the tangent points of the particle velocity to the freeze-out hypersurface in the t-x 
plane. 

For the concave FHS (see Fig. 6.) the freeze-out picture is more complicated because 
there might be a feedback of the particles emitted earlier. In this case the presence of such a 
cut-off term O (p^da^) cannot, of course, exclude reentering particles from the consideration. 
They are necessary for the balance of the energy and momentum and have (in principle) to 
be taken into account. In what follows we shall always mean that only particles which were 
emitted before can reenter the fluid at the concave parts of hypersurface and that, on the 
other hand, particles of the gas with rescattered momenta are not allowed. 

Having this in mind, we can write the distribution function of the gas and each of its 
moments on the FHS denoted by coordinates x* for time t as a sum of two terms 

<j) g (x*, t, p) = (</) g .ermt (x*, t, p) 6 (jfdaj + <f) g . fback (x*, t, p) (-p^da^) , ( 2.12) 

N£ g (x*, *) = J ^ U emU ( x * ? tj p) e + ^ cg fback ( x *, t, p) (-p^da,)) , ( 2.13) 

J p 

(x*, t) = / — V»f U 9 .emit (x*, t, p) (p^daj + 0,. /feacfc (x*, t, p) {-p^da,)) , ( 2.14) 

where the first terms in the right-hand side correspond to the emitted particles, but the 
second ones describe the feedback of the previously emitted particles into the fluid, the vector 
dcr^ being an external normal to the FHS with respect to fluid. We shall also assume that the 
feedback part vanishes if FHS is convex. The notations are as follows: N{? is the 4- vector of 
the current associated with the chemical potential fi c and T^ v is energy-momentum tensor of 
the gas. Throughout this paper we have also adopted that the distribution functions with the 
additional subscript " c" denote the difference of the particles and antiparticles contributions, 
whereas the other ones stand for their sum. 

The above expressions show a main difference between the gas of free particles and the 
usual fluid. One might think of the fluid distribution function cast in the way similar to Eq. 
(2.12) for some hypersurface with a normal vector da^. Formally it is true, but the meaning 
is very different. In the element inside of the fluid the particles with positive (negative) 
product p^do^ are leaving (entering) this element for (from) neighboring element of the 
fluid. But in the gas of free particles there is no neighboring element on one of the sides, 
and only the reentering particles are coming from the earlier emitted gas elements which 
may be quite remote. 

Such a situation should always happen at the time-like boundary of the fluid with the 
vacuum because the traditional hydrodynamics is not applicable at the sharp boundary. 
Hence, at the boundary one has to introduce the "transition matter" that differs from the 
fluid. The main difference between the usual fluid and the "transition matter" is due to the 
cut-off distribution function (2.12) at the boundary with the vacuum. Then one should 
not allow the fluid to have the boundary with vacuum, but should allow the "transition 
matter" to exist between fluid and vacuum on the time-like parts of the boundary. Thus we 
call the emitted particles or the "transition matter" as the gas of free particles in order to 
distinguish it from and to stress its principal difference with the fluid on the time-like FHS . 
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Eqs. (2.12) - (2.14) are the "compromise" between the hydrodynamics and kinetics, 
based on the approximation of a sharp boundary. In sections 4.2 and 4.3 we shall consider 
some justification of this approximation. Hereafter we shall also approximate the distribution 
functions of the emitted particles <fi g .emit by the the equilibrium ones. The latter is consistent 
with the hydrodynamical approach, which implies also the thermodynamical equilibrium. 
Then the adopted assumptions do not seem to be very restrictive because in reality the 
transition region for the time-like FHS has to exist, where, we suppose, there is a fast 
transformation from the fluid into the gas, and due to sufficiently high rate of collisions 
inside the transition region the distribution function of the gas of free particles may be close 
to the equilibrium one. 

The correct formula for one-particle spectra should be applied to the calculation of the 
HBT-correlations from the shock-like freeze-out boundary. First it was done by Sinyukov 
fj~3|| with his one-particle spectrum, and recently the comparison of the Cooper-Frye and 
Sinyukov's formula was made by Heiselberg |14| under very simple assumptions about FHS . 
The same calculation with the correct formula should give new estimations for the observ- 
ables. 

The calculation of the feedback contribution in a closed form is a very interesting task, 
but it is, however, beyond the scope of the present work. We would like only to mention 
some important results with respect to the emission part. In the previous work we gave 
the general expression for the energy-momentum tensor in the RFG and wrote explicitly its 
components for the massless non-interacting gas for the left hemisphere as an example. In 
the subsequent work J| on the similar subject the expressions for the case of massive non- 
interacting Boltzmann gas were presented. Unfortunately, all their formulae are valid for the 
left hemisphere only, and are lacking certain terms in the expressions for energy-momentum 
tensor and the 4-vector of particle number. These results were corrected in part in Ref. 
|J, but in terms of the newly adopted definition for )C n (a, b) functions. Therefore, we would 
like to give in this work a brief technical summary on how to calculate such integrals for the 
time-like FHS which is presented in Appendix A. 

3. Conservation Laws on the Surface between Fluid and Gas 

In the previous section we discussed the derivation of the correct formula for the momen- 
tum spectra on the time-like FHS . From the example of the Sun-like object we saw that 
the Cooper-Frye formula gives wrong results when it is applied for the time-like FHS . This 
property is known for a long time, but the CF formula was used so often in the past because 
it is consistent with the conservation laws of usual relativistic hydrodynamics. However, we 
think it is so because of the mutual cancelation of two mistakes: first, the fact of the decay 
of the small element of the fluid on the time-like FHS has to affect its further evolution (a 
recoil problem, hereafter) and, therefore, has to be taken into account into the equations 
of motion of the fluid, which is ignored in the traditional approach, and, second, for those 
FHS the CF formula of the particle spectra gives incorrect result, but in such a way that 
both mistakes "compensate" each other and hence both are consistent with each other! 

We feel, nevertheless, that it is worth to discuss this problem in more details before the 
formulation of the right approach. 
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3.1. Energy-Momentum Conservation in the Traditional Hydrodynamics and 

the CF formula 

We would like to consider a simple model of freeze-out assuming that the system consists 
of the fluid and the gas of free particles. The transition region between them is not consid- 
ered in this paper. Then fluid evolution is described by the energy-momentum and charge 
conservation equations, or equations of motion in a differential form 

d,Jf(x,t) = 0, (3.1) 
d,N^,t) = 0, (3.2) 

with the evident notations for the energy-momentum tensor T^ v and the charge 4-current 
N"t of the perfect fluid, which are given by the standard formulae 

Tf(x,t) = ( £f +p f )uy f -p f g^ , (3.3) 
Kf{*S) = n c. f u» f , (3.4) 

where Pf, £f and n c ./ are the pressure, energy and charge density of the fluid, respectively, 
and the 4-vector of the hydrodynamic velocity is denoted as = v). 

System of these equations has to be completed by the Additional Conditions such as 

• i) the initial conditions given on the initial hypersurface Sj n (x, t); 

• ii) the equation of state (EoS) of the fluid; 

• iii) the freeze-out criterion, which defines the equation of the FHS denoted as S f r (x, t) . 

We shall assume that the freeze-out criterion follows from some requirement beyond the 
hydrodynamics (for the discussion see, for instance, Refs. [|17| |18|1 and references therein) 
and we shall use just the simplest choice as an example^ to be specific 

T 9 (x,t) = T* = Const => {x*} x 1 * = x\x 2 , x 3 , t) , (3.5) 

i.e., that decay of the gas of free particles happens at the given freeze-out temperature T*, 
and the solution of the left Eq. (3.5) exists. In principle the relevant freeze-out criterion 
should follow from the problem under consideration. For our purpose, however, it is sufficient 
to assume its existence for the gas of free particles and hence it will define the form of the 
FHS . 

We shall denote the latter as x* in the general case, or as x 1 * with the corresponding 
arguments, when it is necessary to distinguish the 4- and 2-dimensions. 

For the sake of simplicity we shall also assume that the FHS is x*, a G and it 

consists of two single- valued parts which are denoted as "left" and "right" ones (in Fig. 5, 
for example, the initial hypersurface S in is AB, and two parts of the FHS one, Ej r , are 
AEF and BGF, respectively). 

2 In the previous paper Q we also considered such criterion as an example only. It is evident that for 
any other choice of the freeze-out criterion , the derived equations can be obtained in the similar way. 
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Let us prove the following Statement 1: Conservation laws in the integral form on the 
FHS follow from the differential ones (3.1) — (3.2) and Additional Conditions if the Cooper- 
Frye formula for the spectrum of the emitted particles is used. 

First we recall the definition of the 4- momentum from the kinetic theory ]TB in terms of 
the momentum distribution function 



P h 



=/ d "V (3 - 6) 



Identifying the gas of free particles with the fluid on the FHS , one can substitute the 
Cooper-Frye formula into the previous equation and thus can obtain the four-momentum of 
the gas of free particles in terms of the fluid energy-momentum tensor 



P u ,fr = J^fp u J *f[ € jr) = J <MT *) > ( 3 - 7 ) 



fr ^fr 



where we interchanged the order of the integrations and substituted then the expression for 
the ideal fluid energy- momentum tensor in terms of its distribution function <j)f. 

Next we integrate the fluid equations of motion (3.1) over the 4- volume of the fluid Vf 
surrounded by the closed hypersurface £/ = Ej n © E/ r with the help of the Gauss theorem 



J d 4 x d^Tf (x, t) = j> da^Tf (x, t) = 0, ( 3.8) 



transforming it into the integral over the closed hypersurface £/. The last step of the proof 
is to rewrite the above integral over closed hypersurface as a sum of two over the initial and 
the final hypersurfaces T> in and £/ r , respectively 

Pf,n = ~ J (*» *) = / *) = P 9-fr , ( 3.9) 

where we have changed the sign in the first equality because da^ is the vector of the ex- 
ternal normal with respect to the fluid. The above equation is just the energy-momentum 
conservation of the fluid in the integral form. 

For the conserved current (3.4) the proof is similar. 

We have to stress that in such a consideration there is no difference between the fluid and 
the gas of free particles . Usually this fact is understood implicitly. In the previous section 
it was shown that the gas of free particles has the cut-off distribution function and, hence, 
differs from the fluid. Thus, the identification of the fluid and gas for the time-like parts of 
the FHS leads to the incorrect freeze-out which, nevertheless, is formally consistent with 
the energy-momentum conservation in the traditional hydrodynamics. Due to the same 
reason, there is no recoil of the particle emission from the time-like parts of the FHS onto 
the fluid evolution. This is seen from the fact that the equation for the FHS (see the 
freeze-out criterion after the identification Tf — T*) does not depend on the expansion of 
the fluid and therefore it can be formulated independently from the equations of motion , 
which should not be the case for the time-like hypersurfaces. Moreover, we know that only 
a part of the particles from the fluid can be emitted from the time-like hypersurfaces and be 
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detected then, and the rest of particles should reenter the fluid. Both of these facts have to 
be taken into account. 

On the other hand, Statement 1 suggests also that the traditional hydrodynamics is not 
consistent with the cut-off formula for the particles emitted from the time-like hypersurface. 
Thus, we are left with the only possibility to use the modified hydrodynamics together with 
the correct formula for particle spectra. 



3.2. Hydrodynamics with Specific Boundary Conditions and the Cut-off 

formula 



Following previous work 0, we write the energy-momentum tensor of the whole system 
consisting of the sum of the fluid and the gas of free particles terms given by Eqs. (3.3) and 
(2.14) , respectively. In order to clarify the presentation, we use a very simple example of the 
FHS with left and right boundaries only as it was introduced in the previous subsection (see 
also Fig. 5 for the 2-dimensional case). Then total energy- momentum tensor and charge 
current read as 



rriflU 

1 tot 



.1* 



.1. . 



X 



N? ot = N£ f e 



i* i 
xi - x 





e 



X 



.1* 



1 1* 
X — XT 



+ 77 {e 



X 



,1* 



•i-- ... 



1 _ 1* 



+ e 
+ e 



.1* 



X 



1* ] 
Xj — X 



} , ( 3.10) 
} , (3-11) 



where the unimportant arguments are dropped for simplicity. 

The meaning of this expression is obvious - it states that the whole system consists of 
the fluid which is " surrounded''^ by the gas of free particles and on the boundary between 
them, x 1 *, there is exchange of energy and momentum. When the hydrodynamical evolution 
is over, there is only a gas of free particles. A similar consideration can be trivially done 
for the conserved currents, but we do not repeat it just to reduce the complications for the 
reader. In what follows it is assumed that the derivatives of the fluid and gas hydrodynamical 
quantities are continuous and finite everywhere in the corresponding domain including the 
FHS. 

The equations of motion of the whole system are the energy-momentum and charge 
conservation laws in the differential form: 



d,TZ = 0, (3.12) 
d,N? ot = 0. (3.13) 

Let us study the boundary conditions first. 

Integrating equations (3.12) over the 4- volume around the vicinity of the FHS , namely 
x 1 G [x 1 * — 5 2 ; x 1 * + 5 2 } and the corresponding finite limits for the other coordinates (see 
parts BDE and BCE, respectively, on Fig. 7 ), and applying the Gauss theorem as it was 
done in previous subsection, one obtains the energy-momentum conservation on the part S f r 
of the FHS 

3 For the time-like parts of the FHS it is true, because there is a spatial separation between fluid and 
gas. But for the space-like ones they are separated in time. 
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3.14) 



-<5 2 



h6* 



where in both sides of the equality do^ is the external normal with respect to the fluid. 

Then in the limit 5 — > for the energy-momentum tensor defined by the Eqs. (3.10) , 
one easily gets the conservation laws for both the energy-momentum and the charge on the 
FHS as follows 



da^Tf" 



or, writing it explicitly, 



d 3 p 



3.15) 
3.16) 

3.17) 



X L * 



/dp / 
p^p u [(pg.emit (p) 6 (p^da^) + (f) g . fback (p) 6 (-p^da p 



da,, 



Po 

d?p „ , ( P P u p f 



X 1 * 



dV M / — p M (0 c . s . emit (p) (p^da^) + (pc.gjback (p) © {-p^da^ 
J Po 

where we have dropped the space-time dependence in the distribution functions above for 
the sake of convenience. 

From the expressions (3.17) and (3.18) it is clear that there are two distinct cases, namely 
(i) when the distribution function of the gas of free particles coincides with that one of the 
fluid (or the EoS is the same for both) and (ii) when they are entirely different. In the 
latter case there is no criterion which equation of state is preferable. This should be defined 
by the physics of the considered task. We believe that the full solution of the problem can 
be given within the kinetic approach only. However, since we adopted the hydrodynamical 
approach which implies also the thermodynamical equilibrium, in what follows it is assumed 
that at least the emission part is described by the equilibrium distribution function. 

Statement 2: If the fluid and the gas have same EoS, then on the space-like parts of 
the FHS their distribution functions are identical, on the convex parts of the time-like 
FHS there exists a freeze-out shock transition, and on the concave ones there is a parametric 
freeze-out shock transition (see below). 

Let us demonstrate this. Indeed, a simple analysis shows that for the case (i) there are 
the following possibilities: 

(i.a) if (p^dan) = 1 for V p, i.e., on the space-like parts of the FHS or on the light 
cone, there is no feedback and then there is only a trivial solution given by the CF -formula 

with the fluid temperature being equal to the freeze-out one. 

Note that the discontinuity through such a hypersurface (or the time-like shocks, as the 
authors of Ref. 



V T 



3.18) 



T 



'g.emit \P) 



19 1 call it) is impossible in this case as it was shown in Gorenstein's work 



0- 
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(i.b) for the convex time-like parts of the FHSthe feedback term vanishes as we dis- 
cussed before (see also Fig. 5) and hence there is a freeze-out shock between the fluid and 
the gas of free particles (actually, the same fluid, but with the freeze-out temperature and 
CO distribution function; for details see later). It is reasonable to call it this way because 
the pressure and the energy density of the gas of free particles are not the usual ones, but 
have extra dependence on the parameters of the FHS . The hydrodynamic parameters of 
the fluid should be found from the conservation laws on the discontinuity 



(PpU P f\ 



da,, 



Po 



P* 1 <t>c.f 



da,, 



da,. 



Po 
Po 



p»p u , 



p p u 



9 (p^da, 



9 (p»da, 



,(3.19) 
. ( 3.20) 



X 1 * 



(i.c) for the concave time-like parts of the FHS the feedback term does not vanish 
(see also Fig. 6) and hence it is a new type of shock between the fluid and the gas of free 
particles which can be called a parametric freeze- out shock with the contribution of the 
feedback particles being a parameter in the general meaning. In this case the equations on 
the discontinuity, i.e. Eqs. (3.17) and the similar one for the conserved charge, should be 
studied. However, it requires the knowledge of the expressions for the feedback particles 
emitted from arbitrary hypersurface which is outside the scope of this paper. In addition, 
the assumption of the instant thermalization of the feedback particles has to be investigated, 
but this is far beyond the hydrodynamical approach. 

Statement 3: If the fluid and the gas have entirely different EoS, then on the space-like 
parts of the FHS time-like shocks |19|| can exist, on the convex parts of the time-like FHS 
there is a freeze-out shock transition, and on the concave ones there is a parametric freeze- out 
shock. 

The proof is similar to the previous consideration. If the fluid and the gas EoS are 
different, [case (ii)], one has the following options: 

(ii.a) if (p^da^) = 1 for V p, i.e., on the space-like parts of the FHS or on the light 
cone, again the feedback contribution is zero, but there is no trivial solution in contrast to 
the case (i.a). And the time-like shocks [|T^] are possible as it was discussed in Ref. pH- 
They are defined by the following conservation laws 



da,. 



da„ 



Po 
Po 



P (p c .f 



T 

f pp u T 



da,, 



da,, 



d*p 
Po 
d 3 p 

Po 



p^p u <p 



P4>c. 



g.emit 



g.emit 




3.21} 



3.22) 



where the equilibrium distribution functions for the gas of free particles 4> g . e mit ( S t*~^) anc ^ 

(pc.g.emit ( E f^) are used. 

More recent results on this subject can be found in Refs. |22fl . However, this kind of 
solutions can probably exist only under very special conditions, namely for the supercooled 
quark-gluon plasma. It is difficult to imagine the reason for the fluid near the freeze-out 
state (when all interactions between particles are very weak) to convert suddenly without 
any cause into the gas of free particles with entirely different EoS . 
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(ii.b) for the convex time-like parts of the FHSthe feedback term again vanishes and 
therefore there is a shock transition between the fluid and the gas of free particles , but now 
fluid and gas being completely different states! Outside the fluid, the gas should be described 
by the CO distribution function. The hydrodynamic parameters of the fluid then are defined 
by the conservation laws on the discontinuity 



da^ 



, ( 3.23) 
. ( 3.24) 



(ii.c) for the concave time-like parts of the FHS the feedback term does not vanish and 
therefore the parametric freeze-out shock introduced above should exist. This kind of the 
shock solution should satisfy the conservation laws in the most general form of Eqs. (3.17) 
and (3.18) . However, the detailed consideration of that solution would lead us too far off 
the main topic of this work. 

Thus, we have given the full analysis of all possible boundary conditions. The con- 
servation laws equations discussed above have to be solved together with the equations of 
motion for the fluid in order to find out the FHS . Next subsection is devoted to the deriva- 
tion of the equations of motion of the fluid alone and their consistency with the boundary 
conditions. 

Note that the conservation laws on the boundary between the fluid and the gas of free 
particles discussed in || (see Eqs. (6), (7) therein) were just postulated and are not related 
to any hydrodynamical evolution of the fluid at all. Moreover, those equations, in the integral 
form as presented in Ref. [|J, cannot be used to solve them together with the hydrodynamical 
equations of the fluid. In this sense those equations are ad hoc ones, and thus of a rather 
restricted use. 

To summarize, we have presented above the full and complete analysis of the possible 
boundary conditions, following the original idea of the paper M. The derived formalism 
allows us not only to find out the new class of the shock transitions, i.e. the parametric 
freeze-out shock, but also to formulate the hydrodynamical approach in the way consistent 
with the emission of the particles from an arbitrary FHS . 



3.3. Consistency Theorem 

Let us study the consistence of the equations of motion for the whole system with the 
boundary conditions derived in the previous subsection. For our present purpose the ex- 
plicit form of the energy-momentum tensor and the conserved 4-current for the gas of free 
particles is not important, but we should remember that those quantities have to satisfy the 
conservation laws in the differential form, i.e., 

= 0, (3.25) 
W£ 9 = (3.26) 
in the domain where the gas of free particles exists. 
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Exploiting this fact, we can rewrite the original conservation laws (3.12) and (3.13) as 
equations of motion of the fluid alone 



x 



1* 



X 



e 



x — X 



1* 



T } = - E * 0# - xl ) n M x (3.27) 



/ —p^p" Ug.emit (P) © (p^/a)) + <Pg.fback (p) © (-p M ^(/a 

Po 



X 



3.28) 



adLr 



f—P^ Uc.g.emit (p) © (p%(/a)) + c . 9 . /feac fc (p) (-?%(/«,))) 



with n^(fa) being the external normal 4- vector with respect to the fluid and being defined 
as follows 

— 1 , a = I , 

n,{fa) = (-dox\*; 1; -c^*; -<9 3 ^*) ( 3.29) 

, +1, a = r, 

i.e., having the positive projection on the X-axis for the right hemisphere and the negative 
projection for the left one (see also Appendix A and Fig. 8.). 

Equations (3.27) and (3.28) look like the hydrodynamical equations of motion with the 
source terms, the first of them describes the loss of the energy-momentum flux through 
the FHSdue to emission of the "frozen" particles and the second one is responsible for 
the gain by the reentering particles at the concave parts of the time-like FHS . It is evi- 
dent because the G-function in the first term takes into account only positive values of the 
product p^n^fa) > 0, while the second one is nonvanishing only for its negative values, i.e., 
P^n^fa) < 0. 

However, it is easy to show that there are no source terms in fact because they disappear 
due to the boundary conditions studied in the previous subsection! Indeed, taking derivatives 
of the remaining 0-functions in the left hand side of Eqs. (3.27) and (3.28) , one obtains 
source-like terms with the fluid energy- momentum tensor and with the conserved 4-current, 
respectively, 
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i* i 

X r - X 



e 



i i* 

X — XT 



W - E 6 { X a - M/a)T u = " E * ( X a ~ >< 



n^(fa) I —P^P U Ug.emit (p) © (p%(/a)) + <Pg.fback (p) © (-£%(/«, 

J Po 



3.30) 



i* _ i 



e 



i i* 

x — x l 



d,K. f - E s 0a* - x 1 ) n,{f a )K f = - £ 5 (x 1 ; - x 1 ) 

agZ,r x a a£Z,r 



e l 

n^ifa) [ — P^ Uc.g.emit (p) © (p P p(j a)) + 4>c.g.fback (p) © (~P M ?V(./a 

J Vn x 



X 



3.31^ 



which automatically cancel the corresponding contributions from the gas of free particles 
due to the general boundary conditions (or conservation laws) on the FHS given by Eqs. 
(3.17) and (3.18) . Then the equations of motion of the fluid acquire the familiar form 






i* 


-x 1 ' 





' x l x 1 * 





1* 


-x 1 ' 





1 1* 

x — x l 



a rplW 



0. 
0. 



3.32) 
3.33) 
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These equations complete the proof of the following Theorem 1: If the gas of free par- 
ticles with the emission part defined by cut-off distribution and with known feedback part 
is described by the equations of motion (3.25) and (3.26) , then the equations of motion 
(3.12) and (3.13) of the whole system consisting of the perfect fluid and the gas of free par- 
ticles split up into the system of the fluid equations of motion (3.32) and (3.33) , and of 
the boundary conditions in the general form of Eqs. (3.17) and (3.18) on the FHS . 

It is necessary to stress that Eqs. (3.32) and (3.33) have only one, but a crucial point of 
difference with the equations of motion of traditional hydrodynamics, although they look 
very similar. Evidently, a solution of the traditional hydrodynamic equations under the given 
Additional Conditions is also a solution of Eqs. (3.32) and (3.33) . Moreover, the solution 
of the latter has to be understood this way. However, the 0-functions in front of the Eqs. 

(3.32) and (3.33) ensure that the fluid exists inside of the FHS (i.e., in the fluid domain) 
defined by the freeze-out criterion for the gas. Thus, the effect of the particle emission is 
implicitly taken into account into the equations of motion of the fluid alone, and is expressed 
explicitly in the boundary conditions between the fluid and the gas. 

Theorem 1 leads also to the fact that inclusion of the source terms which are proportional 
to the 5-function on the FHS into the equations of motion of the fluid will always require 
their vanishing, unless the derivatives of the fluid energy-momentum tensor or the gas one 
contain similar singularities. This is a very important consequence of the Theorem 1 which 
is easy to understand recalling the fact that singularities like 5-function and its derivatives 
(if they appear) are of different order and should be considered independently. Therefore, if 
these singularities enter the same equality, then singularities of the same order will generate 
independent equations. 

Now we are in a position to consider the consistency problem of the energy-momentum 
and the current conservation of the fluid given by Eqs. (3.12) and (3.13) with the emission 
of the "frozen" particles of the gas described by the CO -distribution function like it was 
shown for the usual hydrodynamics and the CF -formula. 

Theorem 2: Energy, momentum and charge of the initial fluid together with the cor- 
responding contribution of the reentering particles from the concave parts of the FHS are 
equal to the corresponding quantities of the emitted from this hypersurface free particles 
with the cut-off momentum distribution function, i.e., those quantities are conserved. 

Proof. First we apply the results of the Statement 1 to the derived equations (3.32) and 

(3.33) of the fluid evolution. Evidently one can use it. Then we obtain an integral form of 
the energy-momentum conservation , for instance, for the fluid alone, which is identical to 



Due to the boundary conditions (3.17) on the FHS which are just conservation laws on the 
boundary between fluid and gas, one obtains the following equality 



Eqs. (3.9) 



P. 



f.in 



J d<Tf(x,*) = J da,T^(x*,t) 





( 3.34) 
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which actually states the energy- momentum conservation in the integral form and also shows 
that the sum of the energy- momentum of the fluid and the particles reentering the fluid during 
its evolution are exactly transformed into the energy- momentum of the emitted particles with 
the correct distribution function! 

A proof for the conserved current can be obtained in a similar way. 



4. Equations of the Freeze-out Hypersurface 

These equations for the time-like parts of the FHS follow from the boundary conditions 
(3.17) and (3.18) . They have to be solved together with the solution of the usual hydrody- 
namical equations of the fluid and allow us to find the FHS as a solution of the system of 
the coupled partial differential equations. It is worth to make a close investigation of these 
equations due to necessity of further applications in both academic and numeric aspects. 

As a good example of applying the derived scheme, we consider an important problem 
of relativistic hydrodynamics - the freeze-out of the simple wave. 



4.1. Freeze-out Calculus: General scheme 



General expressions for the corresponding shock solution on the boundary can be obtained 
from Eqs. (3.17) and (3.18) by the substitution of the formula (3.29) for the normal vector 
in terms of the derivatives of the FHS . In what follows we shall neglect the contribution of 
the feedback particles in order to simplify consideration. Then we shall not distinguish the 
freeze-out and the parametric freeze-out shocks. 

Let us obtain now the full system of all necessary equations. First, we shall suppose that 
the FHS exist and then we derive its equations. For this it is convenient to evaluate the 
boundary conditions in the rest frame of the fluid (hereafter RFF ) before the freeze-out 
shock , where the energy-momentum tensor of the fluid is diagonal. We choose the local 
coordinate system with one of the axes, let it be the X-axis, being parallel to the normal 
3-vector. The latter then is reduced to 



-1 , a = I , 

(/.) = (-fta>; 1; 0; 0) \ (4.1) 

+ 1 , a = r . 



Next, we mention that the velocity of the gas of free particles cannot have nonzero 
projection on any tangential direction to the normal 3-vector in this frame. It can be 
shown directly by manipulation with formulae, but it is evident from the simple reason that 
an oblique shock (see corresponding chapter in Ref. ||23|| ) should have continuous tangent 



velocities on the both sides of shock. Then it follows for the RFF that the gas velocity can 
be parallel or antiparallel to the normal 3-vector only. 

This statement is valid for the RFG as well, and the expression for the normal 4- vector 
in this frame is evidently similar to the Eq. (4.1) . In order to distinguish them, hereafter 
we shall write the corresponding subscript. 
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Boundary conditions (3.17) and (3.18) have a simplest representation in the RFG since 
the moments of the CO distribution function do not look too much complicated in this frame: 

T™{v sG )v sG , ( 4.2) 

T g X0 (v sG )v sG , ( 4.3) 

N° g (v sG )v sG , ( 4.4) 

where the energy-momentum tensor and 4-current of the fluid have standard form of Eqs. 
(3.3) and (3.4) , respectively, with the 3-velocity Vf G . In the above formula the velocity 
v sG = d^x 1 * is the time derivative of the freeze-out hypersurface in the RFG , or velocity 
of the shock in this frame. Note, however, that in contrast to the usual shocks the above 
equations look like the conservation laws in the arbitrary Lorentz frame (not the rest frame of 
the gas!) where the nondiagonal components of the energy-momentum tensor are nonzero. 

Introducing the following notations for the "effective" energy density, pressure and charge 
density of the gas of free particles 

P 9 (v s g) 
n c .g{v sG ) 

one can transform the right-hand side of Eqs. (4.2) - (4.4) to the familiar expressions of 
the relativistic shocks |23|] written in the rest frame of the matter behind the discontinuity. 
In contrast to the usual shock, however, Eqs. (4.2) - (4.4) do not form the closed system 
together with the EoS , but they are dynamical equations for the trajectory of freeze-out. 
Then velocities of fluid and shock in the RFG can be expressed by the standard relations 

(e f -e g (v sG ))( Pf -p g (v sG )) 
(£ f +p g (v sG ))(p f + e g (v sG ))' 
( Pf -p g (v sG )) (e f +p g (v sG )) 
(e f -e g (v sG ))(p f + £ g (v sG ))' 



Tf°(v fG )-T°°(v fG )v sG = T X0 (v sG )~ 
T xx (v fG )-T X0 (v fG )v sG = T xx (v sG )- 
N x (v fG ) - N°Jv fG )v sG = N x (v sG )~ 



= T™(v sG )-T™(v sG )v 7g -, (4.5) 
= T xx (v sG )-T X0 (v sG )v sG , (4.6) 
= N°(v sG ) - N x (v sG )v; G , (4.7) 



v) G = 

VsG = 



Now it is clearly seen that last relation is a transcendental equation for the v sG - velocity of 
the FHS in the RFG . It cannot be solved analytically for an arbitrary EoS . In addition 
we have to transform it to the fluid rest frame in order to complete it with the solution of 
the hydrodynamical equations for the fluid 

V S G -VfG f a tn\ 

V S F = " — • ( 4.10) 

1 - V sG Vf G 

Fortunately, there exist a simple expression for the square of this velocity, namely 

2 (Pf ~ Pg(Vs G )) (Pf + £g(Vs G )) , . ... 

s¥ (e f -e g (v sG ))(e f + p g (v sG )) ' 1 ' > 

which can be easily understood if one recalls that in the theory of relativistic shocks the 
above relation has a meaning of the shock velocity in the rest frame of the initial fluid. 
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Equation for the charge density becomes 



n 



P 9 {v s g) + £f 



~2 / APf + £ f ) 



4.12) 



Evidently, it can be cast in the form of usual Taub adiabate |^4|. Together with equations for 
the shock velocity in RFG and RFF , Eqs. (4.9) and (4.11) respectively, it forms a complete 
system of boundary conditions. 

Let us discuss the boundary conditions and how to solve these equations together with 
the hydrodynamic equations for the fluid. In what follows we shall assume that the solution 
of hydrodynamical equations (3.32) and (3.33) for the fluid is known in the center of mass 
frame (hereafter CM. ) for the whole available space-time volume, and hydrodynamical 
quantities, for instance, T/,/i/ and My C .M. = (1; v /c.m. ) /yl — v /c.m. > are gi ven m eac h 
space-time point XJ C M with v G {0; 1; 2; 3} . Having this solution, one can map it into the 
RFF by the Lorentz transformation 



X F 



1 V /C.M. 



1 - v 2 

1 V /C.M. 



(^C.M. -^C.M. v }c.M. 



^C.M. -^C.M. v }c.M. I ) 



4.13) 

4.14) 



with i G {1; 2; 3} . Then all hydrodynamical quantities are defined in the RFF . 

After this transformation into the RFF , the obtained coordinate system does not nec- 
essarily coincide with the original system that was used for the derivation of the shock-like 
expressions (4.9) , (4.11) and (4.12) on the boundary between fluid and gas. Suppose that in 
the RFF the 3- vector of the shock velocity v s f is described by the standard set of spherical 
coordinates with the angles <p and 0. Exploiting this fact, one can derive differential equation 
for the FHS in RFF as follows. 

First, it is necessary to project the shock velocity v s p on the coordinate axes of the 
RFF . The shock velocity was assumed to be parallel to the X'-axis of the original coordinate 
system. Using the spherical coordinates of the RFF obtained from the CM. frame, one 
can find the projections of the shock velocity v s f on the coordinate follows 



dX 1 



dX° 
dX 2 



dX° 
dX 3 



dX° 



sF 



sF 



sF 



{^f} 1 = v sF sin 9 cos 
{v s f} 2 = v s f sin 9 sin 

{v s f} 3 = V sF COS 9 , 



4.15) 
4.16) 
4.17) 



where superscripts denote respective coordinates. The system obtained above can be used 
in numerical studies of the FHS for the unknown function X*° (X 1 ; X 2 ; X 3 ). 

Next, we reexpress it for the unknown functions X* 1 (X 2 ; X 3 ; X°) as the derivatives of 
the implicit function: 



*i 



OX 
dX° 



sF 



v s f sin 9 cos 



4.1£ 
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dX* 1 
~dX 2 ~ 

dX* 1 



OX 3 



sF 



sF 



tan ' 

• tan 9 cos <j) , 



4.19) 
4.20) 



This is a desired system of equations for X* 1 (X 2 ; X 3 ; X°) which can be employed in numer- 
ical investigations. 

Finally, one can exclude the spherical angles by the simple manipulations and find the 
following equation instead of the system above 



dX* 1 


2 


dX* 1 


2 / g X *l 




2 dX* 1 


dX° 


sF 


1 + dX 3 


sF { + <9X 2 


sF J _ 


~ VsF dX 3 



sF 



( 4.21) 



which has to be solved along with equations (4.9) , (4.11) and (4.12) obtained from the 
boundary conditions, and with the solution of the hydrodynamic equations for the fluid 
which is mapped into the RFF . 

Let us now show how to do this. For that we rewrite all necessary equations in a more 
convenient form using both EoS 



Fnf(T f ,fX f ,fX g ,V sG ) 



Pf 



Pg(v s 



G 



f 



n ig( V sG) (Pf + £g(VsG)) (Pg^sc) + ^(^g)) 



= 



F v sF (Tf,Vf,Vg,V sG ,V sF ) 
F V sG {Tf,Hf,Hg,V sG ) 



J sF 



o , 
o , 



( 4.22) 
( 4.23) 
( 4.24) 



(p f -Pg(vs G )) (Pf + e g (v sG )) 
(e f - e g (v sG )) (e f + p g (v sG )) 

2 (Pf - Pg(VsG)) (g/ + PgjVsG)) 

{£f-£ g {VsG))(Pf + s g {v sG )) 
where the most important arguments are shown. 

The above derivation can be formulated as the following Theorem 3: If solutions of the 
transcendental Eqs. (4.22) - (4.24) exist, then FHS on the time-like boundary between fluid 
and gas is described by the differential equation (4.21) . 

To prove it let us, first, resolve the transcendental equation (4.24) for the mapped solution 
of the hydrodynamic equations into RFF . Suppose it exists and is denoted as 



VsG = v sG (r f (x* l )^ } (x* 1 ),^) 



4.25) 



Substituting it into equation of the charge conservation (4.22) , one finds the relation between 
the chemical potentials of gas and fluid, and fluid temperature. Assume it can be expressed 
as follows 

V 9 = f, g (T f (X* 1 ) , H (X* 1 )) . (4.26) 

Having substitute this solution into Eq. (4.25) and then both equations into expression 
(4.23) , one defines the shock velocity v s p (Tf (X* 1 ) , /j,f (X* 1 )) in terms of the hydrodynami- 
cal solution mapped into the RFF . The last result completes the differential equation (4.21) 
for the time-like parts of the FHS in the RFF . The space-like ones should be obtained in 
accordance with the traditional Cooper-Frye prescription. On the light cones both parts 
should match and this is a requirement to choose the correct root of the transcendental 
equations for the time-like FHS . 
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In the case of zero charge the above consideration simplifies because both chemical po- 
tentials and equation of the charge conservation should be left out. 

Thus, the principal way how to find the time-like FHS is described. Let us now consider 
the freeze-out problem in 1 + 1 dimensions. 



4.2. Freeze-out in 1 + 1 Dimensional Hydrodynamics 

A. Construction of the Solution 

According to the Theorems 1 and 3 we have to conjugate the freeze-out shock with the 
hydrodynamical solution for the fluid. In 1+1 dimensions the original system of equations 
(4.21) - (4.24) for the FHS is greatly simplified. Setting formally dX* 2 ,dX* 3 -> in Eq. 
(4.21) , one finds then in the RFF 



dX* 



dX° 

or, rewriting it in the CM. frame, one gets 

dX* 



= 1^(7) (*•)), (4.27) 

sF 



dX° 



V fC .M. + ^ S F ( 42g) 



sG = ±v 8 f j ~T ( ■ ( 4 -29) 



sC.M. 1 + Vfc.M. V.sF 

by the relativistic addition of the fluid velocity in CM. . 

Dividing Eq. (4.24) by Eq. (4.23) , one obtains the following important expression after 
getting rid of the squares 

,± _^„, e f +p g (v S G 

Pf + £g{VsG) 

Substituting next Eq. (4.29) into Eq. (4.24) , one obtains the second equation for 
velocities 

VsF v f G = ± Pf'PjM , ( 4.30) 

£f - e g {v sG ) 

where we denote possible sign values by the corresponding superscript. 

For the sake of simplicity we shall consider the fluid and gas without charge. Then Eq. 
(4.22) becomes an identity. In what follows we shall impose that the EoS of the fluid is 

l+c| 

Pf = c 2 s afTf Cs = c 2 s Ef, and gas of free particles has the EoS of the ideal gas of massless 
particles p g = ^-T^ = e g /3. 

Such an example is meaningful because intuitively it is clear that once the interaction is 
not important for the gas and its particles are nearly free, then the most natural choice for 
its EoS is the ideal gas one. In order to simplify presentation we consider the massless gas. 
However, this simple choice shows us conceptually most important physical features of the 
freeze-out model. 

The same might be valid for the fluid as well, but we would like to study a more general 
case of the fluid EoS . Then, the case when both EoS are the same is included by the proper 
choice of the speed of sound c s = of the fluid. 
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Using the results of Appendix B, it is easy to obtain the effective energy density and 
pressure of the gas of free particles 

e 9 = -^T*)^r^, (4.31) 

Pg = £ gi T ) Y2 ' ^ ' 

where we employ the results for the right hemisphere. In the further study we shall consider 
only single boundary between the fluid and the gas of free particles , assuming that fluid 
occupied the left hemisphere at the beginning and the gas to be appeared in the right 
hemisphere (in the meaning of the Appendix A). 

It is a remarkable fact that the effective pressure p g for the massless case is exactly the 
same as the pressure in the perpendicular direction to the 3- normal vector (c.f. Eq. (B.4)). 

In this section we shall adopt the freeze-out criterion as it was introduced in Eq. (3.5) 
. We only mention here that one can introduce, for instance, the freeze-out criterion by the 
energy density in the Landau-Lifshitz frame (see Appendix B) or by the effective energy 
density from the equation above. The latter might simplify a consideration essentially. 

Introducing the new variable 

i ^^T >0, (4.33) 

one can rewrite Eqs. (4.29) and (4.30) with the help of the relations (4.31) and (4.32) as 
follows 

„± R + (l-^o) 2 ^^) 

— - ±— , 12 + , 2 , (4.34) 



Rc*-^ 

v2/_ + 



^f^g = ±^ 7-^ , (4.35) 



R + ^ 

sG 



and, excluding the variable R, one finds from Eqs. (4.34) and (4.35) 



( 4.36) 



Using the hyperbolic tangent function one can easily find solutions of Eq. (4.36) . Since 
its solution vf G = 1 is unphysical, we are left only with the following ones 

vf G = =F 2v sF -1±^ - 3 . ( 4.37) 

V sF ' C s 

For the left hemisphere one has to change only the overall sign in the left-hand side of the 
above equation. 
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Then we shall analyze the solution vf G only, because other case can be obtained by the 
substitution v sF — > —v sF . From the inequality \v£ G \ < 1 the available range of the velocity 
v sF follows 

- 1 < v sF < -c 2 s . ( 4.38) 

The maximal value of vf G corresponds to v sF = —c s (see also Fig. 9). 

Inequalities (4.38) show the limiting values of the shock velocity in the RFF which are 
derived by the conservation laws and relativistic causality condition. However, the entropy 
growth condition will give the narrower interval for the allowed values of the velocity v sF 
(see below). 

After some algebra one obtains an expression for the unknown quantity R by substituting 
Eq. (4.38) into Eq. (4.34) or (4.35) 

R+ = (2c 2 + 2v 2 F + v sF (l + c 2 )f(v sF -l) 4 
3 (v sF + Cg) (v 2 F + Cg) 2 

Finally, the result for the unknown fluid temperature on the boundary with the gas reads 
as 



Tt = 



a g (T*) 4 (2c 2 + 2v% + v sF (1 + c 2 s )f (v sF - 1) 



( 4.40) 



3(7/ (v sF + c 2 s ) (v 2 sF + c 2 s ) 2 
The formal solution of the freeze-out problem in 1+1 dimensions follows from the last 
equation: solving it for v sF (Tf(X*)) and integrating Eq. (4.28) with the known hydrody- 
namical solution for the fluid, one finds the desired answer. 

From Eq. (4.40) it is seen that freeze-out criteria T g = T* and Tf = T* are not equivalent 
in general. The only exception is when the FHS is a straight line in (X°,X) plane. This 
occurs in the freeze-out process of the simple wave (see subsection 4.2.C). The criterion 
Tf = T* looks technically simpler because the time derivative of the FHS v sF (T F (X*)) 
is defined by the hydrodynamical solution for the fluid. Then one has to find the gas 
temperature. However, it might be that under the "bad choice" of the freeze-out criterion, 
the gas temperature becomes too low (in the limit v sF (T F (X*)) -> -c 2 s it follows T g -> 0) 
for the applicability of both thermodynamics and hydrodynamics. 

Using the above results, one finds the relative velocity of the fluid in the RFG from Eq. 
(4.8) 



("Jo)* 



+ + ( 4.41) 

where the superscript "+" corresponds to the solution vf G . Causality condition (4.38) plays 
the very same role for the relative velocity of the fluid. 

Dependence of the fluid temperature and the relative velocity of the fluid in the RFG on 
v sF for c 2 s = | is presented in Figs. 10 and 11, respectively. From the Fig. 10 it is 
seen that fluid temperature is always larger than that the gas one. This fact is born in 
a more general statement, namely the energy density of the fluid £/(?/) always exceeds 
the gas one e g (T*) one for c 2 > |. To show its validity we mention only that inequalities 
R+(v sF = -1) = > 1 and > hold for v sF < and c 2 s > §. 



23 



B. Entropy Production in the Freeze out Shock 



Next we would like to study the problem of the thermodynamical stability, or in other 
words the entropy production in the freeze-out shock . This question is of a special interest 
for us because it proves that the freeze-out shock is a new kind of the discontinuity in which 
entropy increases in the rarefaction transition for thermodynamically normal media |25| . 

The fluid entropy flux through the FHS in an arbitrary Lorentz frame is given by 



SfU 



f u f 



da^sf . 



AA2) 
4.43) 



With the help of the previous section the entropy flux in RFG can be written in the general 
form 



SfU^ 



s f 



'(pf-p 9 {v sG )) (e g (v sG ) + Pg (v sG )Y 



4.44) 



(e f -e g (v sG )) (p f + e f ) 

where the normal vector n M acquires the form n p = (—v sG ; 1) in 1 + 1 dimensions. 

The corresponding expression for the gas of free particles has to be found from the cut- 
off distribution function. In the RFG it is evident that the entropy of outgoing particles 
is accounted by the change of the momentum integration volume d 3 p — > Q(p p da p )d 3 p in 
the standard expressions for both classical and quantum statistics. For the Boltzmann 
distribution function 



16| , for example, it yields 
d 3 p 



J ^.jT j [!-]&</>) Q(jp p da») , 



P 
da,, 



4.45) 
4.46) 



The results of Appendix B (see Eq. (B.2) ) then lead to the following expressions of the 
entropy flux of the massless gas 



s v Jv sG ,T) = s g (T*) 



1 - v sG 




dX° n„ s 



4.47) 
4.48) 



where s g (T*) is the entropy density of the gas at the freeze-out temperature. 

One can find then the ratio of the entropy flux on the both sides of the freeze-out shock 



which becomes 



3*4 + 1 



3 v sF + 1 



16a f v* F (v sF 



4.49) 



4.50) 
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for the solution v^ F with c 2 s = |. 

Fig. 12 shows the limiting values of v sF obtained by the thermodynamical stability 
criterion 

> 1 -1 < v sF < -0.479 for a g = a f . ( 4.51) 



,2_1 

'« 3 



The maximal value of the entropy production in the freeze-out shock for the same EoS 
of the fluid is 



max 



i 

4 ( n \ 4 



1.01088 — , ( 4.52) 



i.e., is about one percent for the same number of the internal degrees of freedom of fluid and 
gas. 

Another problem of stability is the mechanical stability (see, for example, Ref. [2~B and 



references therein) of the freeze-out shock. It is of the crucial importance for the freeze-out 
process because it is related to the recoil problem of the freeze-out process on the fluid 
expansion. Usually one might argue that the freeze-out of the small fluid element would 
affect the hydrodynamical solution in the whole future cone of this element. 

However, it is known that in the thermodynamically normal matter the perturbations of 
the small amplitudes propagate with the speed of sound and the rarefaction shocks propagate 
with the subsonic velocity in the local rest frame of the fluid. Thus, one has to conclude 
that a freeze-out shock with the shock velocity in the following range 

- 1 < v+ < - c s ( 4.53) 

is also mechanically stable with respect to the perturbations of the fluid state (c.f. Fig. 12). 
The compression shocks in the fluid we do not consider because on the freeze-out boundary 
there is the energy loss, i.e., it is a rarefaction. Thus, we have proved Statement 4'- For 
the EoS under consideration the perturbations of the fluid state on the time-like FHS are 
slower than the supersonic freeze-out shock and, hence, do not affect the hydrodynamical 
evolution of the fluid inside the FHS . 

Now we see that in contrast to the usual approach the recoil problem is resolved for the 
supersonic velocities of the freeze-out shock in the RFF . This is one of the main goals of 
the suggested freeze-out scheme. However, an investigation of the mechanical stability of 
the freeze-out shock in full requires a separate consideration. 

It is now worthwhile to verify the validity of the adopted approximations by the compar- 
ison of the mean free path A in the fluid and in the gas of free particles . Since A contains a 
cross-section, we rather prefer to study the ratio of the particle density on both sides of the 
freeze-out shock . Such estimation gives the lower limit because of the inequality 



R p 



= p -l<^. ( 4.54) 



Substituting the densities found in the corresponding rest frames (for the gas of free particles 
we use the results obtained in the Eckart frame that are presented in Appendix B), one gets 



' a f (v sF -l) 3 (3v% + l) 
16 a g ( 3 v sF + 1 ) 3 f s 2 F/ 



1 

2 \ 4 



4.55) 
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for c 2 s = |. It is easy to prove that for the same EoS with the same number of degrees of 
freedom the particle density in the fluid is larger than in the gas of free particles, i.e., the 
following inequalities hold (see also Fig. 13) 

> 1 . ( 4.56) 

L s 3 

It would be a problem of internal consistency otherwise, because in the latter case the fluid 
with the larger mean free path should freeze-out before the gas of free particles and the 
whole consideration would become questionable. 

It is yet clear that a very large mean free path in the gas of free particles might lead to the 
reduction of the collision rate and, hence, the usage of the cut-off equilibrium distribution 
function would not be justified. However, such a region of the velocity v s p — > —c 2 s is not 
allowed by the thermodynamical stability condition (c. f. Eq. (4.51) ). Thus, it is shown that 
the considered freeze-out scheme in 1 + 1 dimensions does not have internal contradictions. 
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C. Freeze out of the Simple Wave 

Let us consider the freeze-out of the semi-infinite homogeneous normal matter without 
charge, occupying the left hemisphere in 1+1 dimensions. Then the hydrodynamical solution 
for the fluid is known - it is a simple wave f2~3fl . This simple isentropic flow describes the 



propagation of the perturbations with small amplitudes in the thermodynamically normal 



media |25|. The characteristics of the simple wave and, therefore, its isotherms are just 
straight lines in the space-time variables originating at the initial position of the boundary 
with the vacuum. It is important to remind also that characteristics are the time-like hy- 
persurfaces and, therefore, are of our special interest. A short and clear description of the 
simple wave can be found in the Appendix A of Ref. |EB|. 



Adopting the EoS of the previous two subsections we can apply its results straightfor- 
ward to describe the freeze-out of the simple wave. Making consideration in the RFF , we 
conclude that the shock velocity in this frame, namely v s f, has to be equal to the velocity 
of the simple wave there, i.e., to the velocity of sound c s . Since both waves move to the 
left-hand side in RFF , we write 

v s f = -c s . ( 4.57) 

Hence the freeze-out of the simple wave corresponds to one particular choice of the freeze- 
out shock and it is only necessary to substitute relation (4.57) in all formula of previous 
subsections. 

Thus, substitution of Eq. (4.57) into Eq. (4.37) yields 

vf G = ±i±^ - 3 . ( 4.58) 

The latter equation leads to the restrictions on the fluid EoS . Indeed, requiring the validity 
of the inequality \vf G \ < 1, one gets the lower limit of the velocity of sound in the fluid, i.e., 

c; = 2-V3<c s <l<cf = 2 + V3. ( 4.59) 
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The lower boundary of the velocity of sound corresponds to the extremely soft EoS with 
c~ ~ 0.26795 or with the power of the temperature in the expression for the pressure being 
15. 

The ratio of the energy densities becomes 



l+cj 



s.w. 



'l-Ac s + c 2 s f(l + c s ) 
12 d (1-c) 



( 4.60) 



and the temperature of the fluid in the simple wave reduces to the expression 



s.w. 



Mn 4 (i 



4c s + c 2 s f (1 + c s 



12(7/ C 3 (1-c.) 



el 



4.61' 



In order to illustrate the scale of quantities let us study the case c s = -j= xs 0.57735. 
Then the shock velocity in the RFG is vf G = -j= — 3 « —0.6906, and the ratio R + is 



R+ 



s.w. 



3 



1.5396. The freeze-out temperature of the fluid in the simple wave is 





' 8^ 


s.w. 


_3<T/ \/3_ 



1.1139 



L°7 



4.62) 



If the fluid and the gas of free particles possess the same number of the degrees of 
freedom, then the fluid temperature exceeds the freeze-out one only by 11 percent. Thus, 
no dramatical difference for the temperatures should be expected when both EoS are same. 

The entropy growth in the simple wave corresponds to the maximal possible value 



max Pj 



s.w. 



On 



.2-1 



1.01088 I -± 



4.63) 



The ratio of the mean free paths on the both sides of freeze-out shock satisfies the inequality 



A/ 



s.w. 



> K 



1.655 



s.w. 



( 4.64) 



Thus, the freeze-out scheme of the simple wave is thermodynamically stable and leads to 
the increase of the mean free path in the gas of free particles compared to the fluid. 

An application of the freeze-out of the simple wave in the relativistic nuclear collisions 
should be, of course, done with a more realistic EoS . However, present consideration already 
shows that for the pion rich matter which expands with Landau initial conditions one should 
expect the reduction of the emission volume in comparison with the standard CF freeze-out 
picture. It is so because the same freeze-out temperature of the emitted particles in the 
picture with freeze-out shock and without it corresponds to the different energy densities 
of the fluid in the simple wave. In the freeze-out shock the fluid temperature is larger than 
the freeze-out one, whereas in the traditional approach they are equal. It follows now that 
for the same initial condition the fluid velocity before the freeze-out is smaller for the larger 
energy density, i.e., in the picture with the freeze-out shock. 
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Due to this fact one should expect some reduction of the emission volume for the same 
expansion time for the suggested freeze-out of the simple wave in comparison with the 
traditional Cooper- Frye estimations. This is a clear indication of the "influence effect" of 
the particle emission on the evolution of the fluid in the simple wave. 

Let us now find the spectra of massless particles (bosons to be specific) emitted by the 
simple wave within the cut-off freeze-out scheme and compare it with corresponding result 
obtained by the traditional Cooper-Frye one. Adopting the standard cylindric geometry (X- 
axis is longitudinal) in momentum space, one gets the following expression for the invariant 
spectra in the CM. (after the angular integration) 



dN cut off Et cosh ^ c m ) ( tanh ( yc m ) _ e (tanh (y c .M. ) - v ff ) 

, ( 4.65) 



Ptdptdyc.M. dS ± dt f ^-Hv c .m.-v G )-^ _ ^ 



where E t = \J p\ + m 2 is the transverse energy of the particle, p t is its transverse momentum, 
2/c.m. is the CM. particle rapidity, i] G is the rapidity of the RFG in the CM. frame, and 
v a is a velocity of the freeze-out shock in the CM. frame. For simplicity we consider only 
one internal degree of freedom. 

For the massless particles without charge (/j, c = 0) it is convenient to introduce the 
dimensionless variable p t = Tpr. Then Eq. (4.65) becomes 



T *-3 dN cut off _ p t cosh (yc.M. ) (tanh (y c .M. ) - Vg) 9 (tanh (y c .M. ) - Vg) 
Ptdptdyc.M. dS ± dt " (2n) 2 ( e P* cosh (yc.M. -vg) _ i) 

( 4.66) 

Let us assume that initial fluid has the temperature T in . Then the velocity Vf and the 
temperature Tf of the fluid in the simple wave are related through the expression 

v f (Tf) = tanh (V In . (4.67) 



The freeze-out shock velocity in the CM. frame is given by the formula 

_ V f( T f)- C s 



4.68) 



1 - c 8 v f (T f ) ' 

where the fluid temperature is taken on the freeze-out boundary with the gas and it is 
defined by Eq. (4.40) . 

Substituting v s f = —c s into Eq. (4.41) , one gets the RFG velocity in the CM. frame 

ta„h(, ) = f^f^Z£f) . (4.69, 

1 - Vj G (-C,)Vf (Tf) 

Figs. 14 and 15 show these velocities for the fluid EoS function of the initial 

temperature T in . For such an EoS the freeze-out temperature of the fluid is given by Eq. 
(4.62) . Assuming that fluid and gas have the same number of degrees of freedom, one 
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obtains 





cut- 


-off 


v f 


c 2 = 


1 

3 




cut- 


-off 


fG 


c 2 = 


1 

3 



tanh ( In 



T 



1.1139T* 



-.18835 . 



( 4.70) 
( 4.71) 



The traditional freeze-out scheme based on the Cooper-Frye formula has several impor- 
tant differences. First of all, spectrum does not contain the cut-off ©-function 



T 



*-3 



dN 



CF 



p t cosh (yc.M. ) (tanh (y c . M . ) 



— f 4 72) 

Ptdptdy c .M.dS ± dt (2ir) 2 ( e P* cosh (2/c.M. -m) - \) ' v ' > 

Second, the freeze-out fluid temperature coincides with the gas one, i.e., Tf = T*. Thus, 
the freeze-out in the traditional scheme happens at smaller energy density and at larger 
velocity of the fluid in comparison with the one explained above. And, third, there is no 
difference between the rest frame of fluid and the frame of its decay. Therefore, the fluid 
velocity and the relative velocity of the emitted particles read as 



v f 



J fG 



CF 



j_ = tanh (V3 In 



CF 



= . 



( 4.73) 
( 4.74) 



Now it is easy to see that for the same initial temperature of the fluid the hydrodynamic 
motion of the fluid in the CF freeze-out scheme is more developed (see Figs. 14 and 15), 
i.e., 



CF 



(T in ) > V c 



cut-off 



(Tin) i 



whereas the CM. rapidities are nearly equal 



CF 



(Tin) < 1]G 



cut-off 



(T ir 



4.75) 



4.76) 



First inequality above leads to the important consequences: for the CO freeze-out scheme 
(i) the energy emission per unit time from the FHSis larger, and (ii) the CM. rapidity 
interval with the positive values of the particle spectra is more broad. Inequality (4.76) has 
a negligible effect on the difference of the CO and the CF results. Therefore, one immediately 
gets the following inequality 



dN cut-off 



> 



J j\j Cooper— F 'rye 



Ptdp t dy c .M.dS±dt Ptdp t dy c .M.dS±dt ' 



( 4.77) 



which, evidently, holds for the spectra integrated over rapidity or transversal momentum. 

Due to these reasons the integrated spectra in both schemes are very different. Figs. 
16 - 26 show some typical examples of the spectra integrated over the rapidity and the 
transverse momenta for both schemes. Evidently, the spectra integrated in the small rapidity 
window about yc.M. ~ are extremely different - for the high initial temperatures the 
Cooper-Frye formula gives negative particle numbers everywhere. This feature is hidden 
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if integration is carried for positive values of the rapidity. Nevertheless, the quantitative 
difference remains. The spectra integrated over the transverse momenta have, basically, 
similar features. However, these simple examples of the spectra show that some previous 
results based on the Cooper-Frye formula may be considerably revised while the correct 
freeze-out scheme is used. 

An application of Eq. (4.77) to describe the finite systems requires the knowledge of the 
hydrodynamic evolution of the system. Nevertheless, it is clear that for the cut-off freeze- 



out scheme the fluid volume is reduced on the magnitude of about ( S± F (t em ) v£ F — Sl° (t em ) J t, 



in comparison with the CF one (S±(t em ) is the area of the freeze-out 3-surface at the emis- 
sion time t em ). The maximal difference of the velocities is about v^ F — « 0.18 and it 
corresponds to the initial temperature Tj n « 1.5T*. Such a reduction of the emitting source 
size might be in a better agreement with the results of the HBT analysis of the relativistic 
nuclear collisions. 

In the considered example we gave the solution of the simple wave freeze-out problem. 
Due to the fact that characteristics and isotherms of the fluid in the simple wave are just 
straight lines in any Lorentz frame, one can show that the ratio of the partial derivatives of 
the temperatures on both sides of the shock is the same, i.e., = in any frame. This 

is the simplest example against the statement of Ref. || that such a condition is unphysical. 
As was shown in this section, it is unnecessarily strong for the derivation of the equation for 
the FHSas it was supposed in the paper 0, but it appears in particular applications like 
the freeze-out of the simple wave. 
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5. Conclusions 



In this work the solution of the freeze-out problem in relativistic hydrodynamics is pre- 
sented within a zero width FHS approximation. We analyze the difference between the 
CF procedure and its CO generalization for time-like FHS . It is shown that a reformulation 
of the traditional hydrodynamic scheme is necessary in order to include particle emission 
from the time-like parts of the FHS . The modified self-consistent hydrodynamics with the 
specific boundary conditions is formulated and different types of shock-like freeze-out are 
studied. 

The correct boundary conditions enable us to derive the equations of motion of the fluid 
alone, which do not contain any source term. We have also proved the energy-momentum 
and charge conservation in the integral form for hydrodynamics with the specific boundary 
conditions if the CO distribution function is used. 

We analyze the equations for the time-like parts of the FHS and show how to solve them 
together with the equations of motion of the fluid. A complete analysis of the freeze-out 
problem for the time-like parts of the FHS in 1+1 dimensions is presented. The entropy 
growth and mechanical stability conditions as well as the "recoil problem" for the shock- 
like freeze-out are studied. As an application of the general scheme, the freeze-out of the 
simple wave is considered, and analytical solution of this problem for the massless gas of free 
particles is given. The spectra of the emitted particles are calculated and compared to those 
obtained in the CF procedure. 
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Appendix A 



The moments of the cut-off distribution function in 4-dimensions can be easily calculated 
in the RFG . This is a convenient frame for such a calculation. 



A.l. Spherical Coordinates in Momentum Space 

Let us assume for simplicity that in 4-dimensions the fluid occupies the left semi-infinite 
volume (see Fig. 8, case a) and its FHS has the following external normal vector with respect 
to the fluid Q 

da* = (-v a ; 1;0;0) dtdydz, (A.l) 

with v a = ^jft having the same meaning as in Eq. (2.11) . From its explicit form it is 
clear that the above vector is always pointing outside the fluid (or to the right hand side of 
Fig. 8.) which we indicated by the Roman superscript in capital. 

Then one can calculate, for example, the energy-momentum tensor in the RFG using 
spherical coordinates in momentum space. Evidently one can explicitly perform angular 
integration if the distribution function depends on the energy and not momentum. This is 
true for any equilibrium distribution function of the gas of free particles in the RFG because 
it is the rest frame of the original, non-cut, distribution function. However, the presence of 
the cut-off B-function makes it more complicated in general, because integrations over the 
energy and spherical angles are not decoupled. The only exception is the case of massless 
particles when integrations over the energy and spherical angles are independent as we saw 
it during the discussion of the emission spectrum for the Sun like object. Nevertheless, one 
can try to find the closed form for such integrals. 

Now if we choose the spherical angle 8 P between X-axis and the vector of 3-momentum 
p, then the integration over this angle is reduced to simple quadratures as follows 



I^{v a ,p) = jd9 p sin(e p )F{cos(0 p ))e(^cos(e p )-v a )= jdxF(x)e(^x-v^ , 

-1 

( A.2) 

where we use the evident notations p = |p|, E — \fp 2 + m 2 . 

The geometrical meaning of this result is clear from the case a) of Fig. 8. The O (p^da^j- 
function cuts off the projections of the 3- velocity on the X-axis direction from the lower limit. 

Similarly, for the fluid occupying the right semi-infinite volume (see Fig. 8, case b)) with 
the external normal vector to the fluid 

da% = (v ff ;-l;0;0) dtdydz ( A.3) 

one has to calculate similar integrals 

4 Note that this is not a very restrictive choice of the normal vector because the general case can be 
reduced to the considered one by the proper rotation of the coordinate system in each space point. 
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7T +1 

I^(v a ,p) = Jd9 p sm(9 p )F(cos(9 p ))e(v a -^cos(9 p )) = J dx F(x)Q(v a - |x) . ( A.4) 



o -i 
These integrals are connected by simple relations, for example, 

I${v a ,p) = I$(+l,p) - I$(v a ,p) , 



(A.5) 



and therefore knowing one of them is enough to find out the other one, and, hence, we shall 
give it just for the left case and then write it for the right one without detailed calculation. 



A. 2. Angular Integration 

In order to find the angular integrals let us first change the variables 

+i 



= I dx F(x)e(v a -L x ) = - [ dr F(-r)e(v a - r) . 
J hi p J p 



,E 
p 



A.6) 



Then it is evident that the integral does not vanish only in two cases, i.e., when v a > and 
— -| < v a < -|, which can be written in terms of 6-functions as follows 



£ J dr F(2t) , v a > 

p j \ p i i o 



+v a 



f / drF(fr), -|<^<| 



p_ 

E 



v a < -f 



A. 7) 



ok- | 



+1 



x 



e(v a ) j dx F(x) + e^-v^e(v a + ^j j dx f( 

Now the integral is decoupled into two parts, the first of them looking like a space-like 
contribution because of the integration limits, and the second one reminds the time-like 
one! One can call it this way since in the limit of space-like FHS (v a > 1) all values of the 
spherical angle 6 P contribute and, hence, only the first term survives. 

However, the meaning of the first term above is deeper: if FHS moves inside the fluid 
with the velocity v a and one considers the values of the particle velocity with the modulus 
-| < then particles with any projection of the velocity ^ are emitted from the fluid, and 
this fact is taken into account by the limits of integration from —1 to +1. For the space-like 
FHS it is always so and the space-like term in Eq. (A. 7) does not vanish. The second term 
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in Eq. (A. 7) , evidently, describes the contribution of the velocity projection on the X-axis 
which exceeds the velocity of the FHS . 

The same is valid for the right hemisphere as one can see from the expression below that 
was obtained with the help of the identity (A. 5) 



+1 +i 



Ip(v a ,p) = (-v a - |) (-v a ) J dx F(x) + (| - v a ) (v a + |) J dx F(x) , 



A.8) 



where we used the fact that for positive -| > the following product of ©-functions always 

vanishes: {-v a - §) (v a - §) = 0. 

Note that similar results hold for the case of the cylindrical geometry, with some simple 



modifications. 



A. 3. Momentum Integration 

Since our next step is an integration over the modulus of momentum or energy, it is 
necessary to rewrite the above results in terms of energy. After some algebra, one gets the 
following result for the angular integrals 



+1 J v « 
l£(v a:P ) = (E a - E)Q (v a ) J dr F(t) + Q (E - E a ) J dr F(r) , 

-l -l 

+i +i 
If(v*, P) = {Ea -E)Q (-v a ) J dr F(r) + (E - E a ) J dr F(r) , 



A.9) 



with E a = y=j f° r v l < 1 an d E a = oo for v 2 a > I. Note that (±v a ) in front of the 

space-like integral above is of a crucial importance when one uses integration over the energy 
on the time-like FHS. For the space-like FHS (i.e., v a > 1) it is trivial because the other 
©-function ensures that the correct value of v a is taken into account. 

In the derivation of Eqs. (A.9) one has to use the following two relations for ©-functions 
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E 



P 



E 



P 



E 



P 



E 



and 



© (| - V^j (V a ) = (E -E'er) © {V a ) , 

from which all other necessary relations can be trivially found. 



© (| - Va) © M + © (v a + |) (~V a ) , ( A. 10) 

(All) 
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Usually one has to integrate such expressions over momentum or energy, i.e., to calculate 
the integrals of the following type: 

oo oo 

l£ F {v a ) = J dpp 2 G(p,E)I$(v a ,p) = J dEEp(E)G(p(E),E)I$(v a ,p(E)), ( A.12) 

m 

with A e {L, R} and p(E) = ^E 2 - m 2 . Using Eqs. (A.9) , one obtains the relations for 
momentum integrals in a closed form, 



(v a ) = 6 (v a ) J dEEp{E)G{p{E),E) J dr F{r) 

m 

+ J dE Ep(E) G (p(E), E) j aJm dr F(r) , ( A.13) 

Ea +1 

IgfM = Q(-v a ) I dEEp(E)G(p(E),E) dr F(r) 



m 

00 



+ ( dEEp{E)G{p{E),E) f dr F{t) . ( A. 14) 

It is evident that a sum of Eqs. (A.13) and (A. 14) gives just the usual result without the 
cut-off 6-function 



00 

4W + 4W = 4(+1)=4H) = JdEEp(E)G(p(E),E) j + drF{r), ( A.15) 

m 

when integrations over spherical angle and energy are not connected. Another example 
of such a simplification is the case of massless particles, when the space-like contributions 
vanish, as it immediately follows from Eqs. (A.13) and (A. 14) 



l GF 



m=0 



00 

= J dEE 2 G (E, E) p dr F(t) , ( A. 16) 



IgfM 





00 



^ ■ / dEE 2 G(E,E) f +1 dr F(t). { A.17) 



m=0 



All these relations can be greatly simplified for even or odd functions F(cos(6 p )) of 
the cosine of the spherical angle 6 P which is usually the case for most applications. After 
trivial manipulations, one can easily prove the validity of the following relations between the 
integrals over left and right hemispheres 

f j gf(-^), even F(t) , 
I§fM = ( A.18) 

I -Igf{-v*) , odd F(r) , 
which is one of the most useful and powerful relation. 
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An important consequence of this statement is that the terms proportional to (±f CT ) do 
not vanish for even functions F(t), and those terms are not present in the final expressions 
for odd functions F(t). This is clearly seen from Eqs. (A. 13) and (A. 14) . 

Using the results above, it is easy to find out expressions for the moments of the distri- 
bution function of the gas of free particles 



A. 4. Moments of the Cut-off Distribution Function 

As a good example to apply the results of the previous subsection, let us consider a 
relativistic massive Boltzmann gas with the distribution function 



E(p)-jt e (x,t) 

(x,t,p) = C g e T(x, t , 
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with the constant C„ 



where d g is the degeneracy number of the gas. Then the 



y 9 (2tt) 3 ' ,ViiVji ^ 

components of the particle flow number of the emission part are 
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for the left hemisphere, and 



Af L K,T,/x c ) = 



-p T 

^2 c 



A.20) 

A.21) 
A.22) 



N g,R( V °> T > 



2TrC n T 3 e^l<S)(-v n 



. m m\ fm E a 
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for the right one. Here we used the following notation for the KVfunction 

2"- 1 (n - 1)! " 



/C n (a; 6) = 



(2n - 2) 



J dxx (x 2 — 



2 e~ x 
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which differs from the notation of the paper and the book [|16[] by the dimensionless factor 
a 11 . 

With the help of the following identity 
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K 2 (0;fW 2 + 2 | + ff)e-* 



A.27) 



we complete our analysis of the particle flow number. 

One can see that in the expression for N gL {y ai T, /x c ) of Ref. f| the factors with Q(v a ) 
and Q(—v a ) are missing, as well as in their expressions for the diagonal components of the 
energy-momentum tensor of the gas of free particles . 

The results for the right hemisphere we give for the first time. 

We also would like to emphasize the fact that the spatial projections of the particle flow 
4-vector on other than the X-axis vanish in the RFG due to convenient choice of the nor- 
mal vector. Same property one sees when evaluating components of the energy-momentum 
tensor . 

Expressions for the nonzero components of the energy-momentum tensor of the emitted 
particles then read as 
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for the left hemisphere. In deriving equations above we used several useful identities 
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j dxx 2 (x 2 -a 2 f e~ x = b 2 {b 2 -a 2 f e- b + 3}C2 (a; b)+ a 2 Id (a; b) , ( A.34) 

b 

°° 3 3 

Jdx (x 2 -a 2 y e - x = (b 2 -a 2 )~ 2 e~ b + ZK 2 (a-b) ; ( A.35) 

b 

The identities similar to these ones were used in Ref. [Q, but the first terms in the r.h.s which 
appear after integration by parts were missing. For the further evaluation of the radicals 

above it is also convenient to use the identity (7^ — ^\ 2 = \v a \^-. 

The formulae for the right hemisphere can be obtained in the similar way by a straight- 
forward integration or with the help of Eq. (A. 18) 
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These expressions can be used for the approximate description of nucleons and baryon res- 
onances in the numerical studies, with the only exception for the very low temperatures 
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where the Fermi statistics is essential. They also show us that even in the simplest case of 
the time-like FHSthe derived expressions are more involved then in case of space-like one. 
The limit of the massless gas follows straightforwardly from the equations above. 



Appendix B 

In case of massless particles the formulae of the previous subsection become simpler. 
Note that for high enough temperatures, namely T > m n , m n is a pion mass, these results 
can be also used as a good approximation for the pion gas. 

B.l. Moments of the Cut-off Distribution Function for the Massless Gas 

Employing Eqs. (A. 16) - (A. 18) one obtains the following results for the particle flow 
number for the left and right hemispheres in the RFG 



n g (T, \i c 



1+Vg 
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l — v„ 



-, 0; 



-; 0;0 
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B.2) 



2 ' 4 

where n g (T, /i c ) is usual particle density of the gas without cut-off with the temperature T 
and chemical potential fi. Similarly, the emission part of the energy-momentum tensor of 
the particles moving antiparallel to the X-axis is given by the formula 



2 4 
vl-1 1+vl 
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(1+^)2(2-^) 
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where e g (T, /i c ) is the usual energy density of the gas with the temperature T and chemical 
potential /i c , obtained for non-cut distribution function. It is remarkable that the above 
results do not depend on whether the statistics is Boltzmann, Fermi-Dirac or Bose-Einstein 
one! The corresponding formula for the case a) of Fig. 8. (or for the particles moving parallel 
to the X-axis) can be obtained similarly or employing Eqs. (A. 5) , (A. 16) - (A. 18) 



T£ t (v a ,T, f i c )=e g (T,/i c 
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(B.4) 



B.2. Eckart and Landau Lifshitz Velocities for the Massless Gas 
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This question was briefly discussed in Ref. 0. We would like to discuss it in more details 
for the gas of massless particles because on this simplest example one can analytically show 
that the value of the Eckart velocity always exceeds the Landau-Lifshitz one. 

First we find out the value of the velocity in the Eckart meaning. Instead of calculating 
it directly from the definition which can be found in text-books (see, for example, Ref. fL6|), 
we remind the reader that in the Eckart frame the particle flow number has only the time 
component nonvanishing. Then, making the proper Lorentz transformation, we immediately 
obtain the expression for the Eckart velocity in the considered case 

N^( VtT ,T,fi c ) = 1 (N^ A (v a ,T,fi c )-v EA N giA (v a ,T,fi c ))=0 ( B.5) 

=> V E)A = -ft— r = I Va \l a _ p B - 6 

Note that such a transformation and, therefore, the first part of Eq. (B.6) are valid for 
massive gas as long as the normal vector to the FHS is chosen in the same way. 

In accordance with the Eckart definition, the proper particle number density reads then 

as 



l+Vt 



■y/(l+V a ) (3-tv), A 



nE,A = n g (T,fi c ) \ ( B.7) 

1^/(1-^) (3 + v a ), A = R, 

which is, evidently, reduced in comparison with the value of n g (T, /i c ) . 

Now we can find the Landau-Lifshitz velocity by the similar Lorentz transformation 
for the non-diagonal component of the symmetric energy-momentum tensor (dropping their 
arguments for simplicity) 

V T, Mc ) = 7 2 ((l + v\ LA ) - v LL , A T™ - v LLA T*j) = ( B.8) 

T^OO I rpXX I ( rpOQ T rpXX r ) rpQX\ I ' rpQQ T rpXX , OT^OXA 

1 g,A + 1 aA - \/ ^ 9 ,A + 1 g,A ~ A1 g,AJ V g,A + 1 g,A + A1 g,AJ , 
=> V LL,A = 1 ' (R9) 

where we have chosen the physical solution of the Landau-Lifshitz velocity by the condition 
vll,a -> when T° g x A -> 0. 

Applying the general result of Eq. (B.9) to the massless particle gas, one obtains more 
involved expressions for the Landau-Lifshitz velocity 
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3(1-d ct ) 
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A — L. 



A = R. 



B.10) 



To compare it with the Eckart velocity, let us find the limit v c — > +1 [v a — > — 1) for the 
left (right) hemisphere. Then it is convenient to rewrite the previous expression in terms of 
the difference A 



40 



vlla 



(2 - A) V4 + 2A + A 2 - (4 - A + A 2 
3A 



A 



-1, A = l + v a , 
Taking the limit A — > 0, one trivially gets the following result 
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Thus, direct comparison with the expression (B.6) indicates that the modulus of the Eckart 
velocity exceeds the Landau-Lifshitz one if the time-like FHS is close to the corresponding 
light cone, i.e., if < A 2 < 1 



v LLA {&) < v e ,a(&) 



( B.13) 



It can be shown that inequality (B.13) holds always for the physical values of the velocity 
v% < 1 for the time-like FHS or for any physical values of A. Indeed, considering only the 
case of right hemisphere, one can find all roots of the equation 



vll,r (A) = v e ,r (A) 



4 - A + A 2 - (2 - A) V4 + 2A + A 2 = -A 2 ( B.14) 
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which are then A = {0; 2}, and only the zero root is physical. Therefore the Eckart and 
Landau-Lifshitz velocities coincide only for the case A = 0, hence inequality (B.13) cannot 
be broken. 



B.3. Anisotropy of the Energy-Momentum Tensor of the Massless Gas 

Results of the previous subsection enable us to study the spatial structure of the energy- 
momentum tensor of the massless particles. The spatial anisotropy of the Cooper- Frye spectra 
was originally found for the space-like FHS ]12[. Here we would like to discuss the spatial 
anisotropy of the energy- momentum tensor on the time-like FHS . The nontrivial struc- 
ture of the gas energy-momentum tensor on the time-like FHS already indicates that such 
anisotropy might exist. Now we can prove it rigorously for the gas of massless particles. For 
this purpose we diagonalize the energy- momentum tensor by the Lorentz transformation of 
Eq. (B.8) , first, calculate the T*Jf component, next, and compare it with the T^Jf one. 

After the diagonalization of the energy-momentum tensor the T** and Tg° A components 
read as 

T g x f(v a ,T,» c ) = 7 2 W(?/-2 W 5+^5). ( B.15) 

T™a(v*, T, M c ) = 7 2 (vll,a) (T 9 °° a - 2v LLA TfX + v\ LA Tfj) , ( B.16) 
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with the Landau-Lifshitz velocity vll,a- It is sufficient to consider just the case of the right 
hemisphere, i.e., A = R. Then in terms of the difference A, introduced in Eq. (B.ll) , the 
energy-momentum tensor in the RFG becomes 
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Employing this convenient form of the energy-momentum tensor and the expression 
(B.ll) for the Landau-Lifshitz velocity, one can find the XX component in the Landau- 
Lifshitz rest frame where the energy-momentum tensor is diagonal. After some algebra one 
gets the following result for Eqs. (B.15) and (B.16) 



(Vl + f + f + l~f) 
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It is only important to mention that in deriving the equations above it is convenient 

A I A 2 
2 "T" 4 



to use the polynomial S(A) = 1 + ^ + in terms of which all expressions are greatly 



simplified, for example, 



1 ± Vll r 



(A), 



^(A) + (l-f) x a 
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Then a direct comparison with the F K component of the energy-momentum tensor from 
Eq. (B.17) indicates that for the physical values of A, i.e., < A < 2, in the Landau-Lifshitz 
rest frame the pressure in X direction is always smaller than in Y or Z one 



Evidently this result is valid for the left hemisphere as well. 



( B.21) 
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Fig. 1. Decay of the element Ax = AB during the time At = BD in the reference frame 
of the gas. Particles cross the time-like hypersurface AD. Trajectories of the particles that 
have negative momenta and are not reflected from the wall BD in the meaning of Ref. 
are shown by arrows. 
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Fig. 2. Decay of the element Ax = AB in the spirit of Ref. Trajectories of the particles 
with negative momentum which are reflected from the wall BD are shown by the solid lines 
with arrows. However, in reality there is no wall, and particles with positive momentum, 
which are emitted from the element AB, are going back to the fluid and their trajectories 
"behind" the wall BD are shown by the dashed lines with arrows. Reflection from the wall 
is equivalent to the consideration of the particles coming from the element BC = AB which 
is located "behind" the wall (trajectories of such particles are shown by the corresponding 
dashed lines). 
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Fig. 3. Decay of the element Ax = AB during the time At = BD. Trajectories of the 
particles that have negative momenta and are reflected from the wall BD in the meaning of 
Ref. |T3 are shown by the dashed lines with arrows. In the real situation those particles are 



going back to the fluid. The correct trajectories of the free streaming particles are shown by 
the solid lines with arrows. Thus one can see that in addition to the contribution from the 
element BC = AB that is taken by the reflection, the contribution from the element CH 
should be taken into account for the momenta — 1 < ^ < ~^r>- ft is only necessary for the 
time-like freeze-out hypersurfaces. 
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Fig. 4. Decay of the element Ax = AB during the time At = BD. Trajectories of the 
particles with positive momenta are shown by the lines with arrows. Such particles will 
cross the decay hypersurface AD only if they are emitted from the element AC = Ax — ^At 

or, equivalently, if the following inequality for the velocities 4f > ^ holds. 



48 




Fig. 5. Decay of the fluid into the gas of free particles on the convex freeze-out hypersurface 
AEFGB. Trajectories of the particles are indicated by the lines with arrows. Dashed line 
represents the detector's world line. Lines AC and BD denote the light cones. For the given 
particle velocity v = ^ the integration limits in coordinate space, points E and G, are found 
from the condition p p da p = 0. The geometrical meaning of points E and G is evident from the 
construction: they are tangent points of the particle velocity to the freeze-out hypersurface 
in t — x plane. In contrast to the cut-off formula for invariant spectra which ensures such 
limits automatically, the Cooper-Frye one takes into account the particles emitted from all 
points of the freeze-out hypersurface AEFGB. 
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Fig. 6. Decay of the fluid into the gas of free particles on the concave freeze-out hypersurface 
ABCEFG. Notations correspond to the previous figure. In this case, however, one has to 
take into account the particles feeding back to fluid on the part ABC of the freeze-out 
hypersurface. These particles were emitted earlier and do not appear from rescattering. The 
latter is forbidden by the assumption that particle spectra are "frozen" once they belong to 
the gas of free particles. Hence, particle trajectories, like the one shown by the dotted line, 
are not allowed in the freeze-out picture because those particles appear from nothing. 
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Fig. 7. Integration contour BCED showing the derivation of the boundary conditions 
between the fluid and the gas of free particles on the FHS ABEFG. 
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Fig. 8. Decay of the fluid into the gas of free particles on the freeze-out hypersurfaces 
ABE and EFG is shown for the right hemisphere (Case a)) and for the left one (Case b)). 
Emission trajectories with possible momenta from point B are indicated by the line with 
arrows. They are restricted by the tanged line BC to the hypersurface and light cone BD. 
For the point F the construction is symmetric with respect to reflection. 
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Fig. 9. Dependence of the shock velocity v s g in the RFG upon the the shock velocity v s f 
in the RFFfor three different fluid EoS . The values of the fluid velocity of sound are as 
follows: c s = -y= (solid line), c s = -4= (dashed line) and c s = 4= (dotted line). 
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Fig. 10. Ratio of the energy densities on the both sides of the freeze-out shock as a function 
of the shock velocity v sF in the RFFfor three different fluid EoS . The legend corresponds 
to Fig. 9. 
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Fig. 12. Ratio of the entropies on the both sides of the freeze-out shock as a function of 
the shock velocity v sF in the RFFfor the fluid velocity of sound c s = 
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Fig. 14. Fluid velocity of the freeze-out of the simple wave as a function of the initial temper- 
ature of the fluid T in . The dashed line corresponds to the Cooper- Frye freeze-out scheme, 
and the solid one corresponds to the cut-off freeze-out scheme. 
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Fig. 15. Velocity of the freeze-out boundary v c in the CM. frame as a function of the 
initial temperature of the fluid T in . The dashed line corresponds to the Cooper- Frye freeze- 
out scheme, and the solid one corresponds to the cut-off freeze-out scheme. The dotted 
line represents the velocity of the RFG in the CM. frame. 
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Fig. 16. Momentum distribution of the Cooper- Frye freeze-out scheme integrated over the 
CM. rapidity yc.M. £ [—1; 1 ]• The spectra are found for the three values of the initial 
temperature of the fluid: T in = LIT* (dotted line), T in = 1.5T* (solid line), and T in = 1.9T* 
(dashed line). 
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Fig. 17. Momentum distribution of the cut-off freeze-out scheme integrated over the CM. 
rapidity yc.M. £ [ — 1; 1 ]■ The spectra are found for the three values of the initial tempera- 
ture of the fluid: T in = LIT* (dotted line), T in = 1.5T* (solid line), and T in = 1.9T* (dashed 
line) . 
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Fig. 18. Momentum distribution of the Cooper-Frye freeze-out scheme integrated over the 
CM. rapidity yc.M. £ [—2; 2 ]. The spectra are found for the three values of the initial 
temperature of the fluid: T in = LIT* (dotted line), T in = 1.5T* (solid line), and T in = 1.9T* 
(dashed line). 
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Fig. 19. Momentum distribution of the cut-off freeze-out scheme integrated over the CM. 
rapidity yc.M. £ [ — 2; 2 ]. The spectra are found for the three values of the initial tempera- 
ture of the fluid: T in = LIT* (dotted line), T in = 1.5T* (solid line), and T in = 1.9T* (dashed 
line) . 
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Fig. 20. Comparison of the momentum distribution of the Cooper- Frye (dashed line) and 
the cut-off (solid line) freeze-out schemes integrated over the CM. rapidity yc.M. £ 
[-2; 2 ]. Initial temperature of the fluid is T in = 1.5T*. 
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Fig. 21. Comparison of the momentum distribution of the Cooper- Frye (dashed line) and 
the cut-off (solid line) freeze-out schemes integrated over the CM. rapidity yc.M. £ 
[-2; 2 ]. Initial temperature of the fluid is T in = 1.9T*. 
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Fig. 22. Comparison of the momentum distribution of the Cooper- Frye 
the cut-off (solid line) freeze-out schemes integrated over the CM. 
[— 1; 6 ]. Initial temperature of the fluid is T in = 1.5T*. 
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Fig. 23. Momentum distribution of the Cooper-Frye freeze-out scheme integrated over 
the CM. transverse momentum p t G [0.01; 6 ]T*. Initial temperatures of the fluid are: 
T in = 1.17* (dotted line), T in = 1.5T* (solid line), and T in = 1.9T* (dashed line). 
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Fig. 24. Momentum distribution of the cut-off freeze-out scheme integrated over the CM. 
transverse momentum p t G [O.Of; 6 ]T*. Initial temperatures of the fluid are: T in = LIT* 
(dotted line), T in = 1.5T* (solid line), and T in = 1.9T* (dashed line). 
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Fig. 25. Comparison of the Cooper-Frye momentum distribution (dashed line) and the 
cut-off one (solid line) integrated over the CM. transverse momentum p t e [0.01; 6 ]T*. 
Initial temperature of the fluid is T in = 1.5T*. 



69 




-0.01 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
-10 12 3 



Fig. 26. Comparison of the Cooper-Frye momentum distribution (dashed line) and the 
cut-off one (solid line) integrated over the CM. transverse momentum p t e [0.01; 6 ]T*. 
Initial temperature of the fluid is T in = 1.9T*. 
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